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1. INTRODUCTION 



Consider the matrix Schrodinger equation on the half hne 



+ V{x) = k^ij, 



X e (0, +00), 



(1.1) 



where the prime denotes the derivative with respect to the spatial coordinate x and the po- 
tential V is a nxn matrix-valued function belonging to class L];(R+) with R+ := (0, +00). 
Note that V e means that each entry of the matrix V is Lebesgue measurable on 

the interval I and 



where ||y(a;)|| denotes a matrix norm. Clearly, a matrix-valued function belongs L];(R+) 
if and only if each entry of that matrix belongs to L\{Il'^). Note that V is not assumed to 
be real valued, but we impose the condition that it is selfadjoint, i.e. 



where the dagger denotes the adjoint (complex conjugate and matrix transpose). Without 
loss of any generality we can view the wavefunction i(;{k,x) appearing in (1.1) either as a 
vector-valued function with n components or as an n x p matrix- valued function for some 
p with 1 < p < n. 

We are interested in studying (1.1) with a selfadjoint potential V in L];(R+) under 
the most general selfadjoint boundary condition at a; = 0. This is the generalization of 
the scalar version (with n = 1) of the corresponding problem, where the most general 
selfadjoint boundary condition at x = can be stated as [5,8,16,30] 



where the parameter 9 takes values in the interval (0,7r]. The special choice 9 = tt corre- 
sponds to the Dirichlet boundary condition, and the choice 9 = 7r/2 corresponds to the 
Neumann boundary condition. 

A formulation of the most general selfadjoint boundary condition at a; = for (1.1) 
was stated in [23,24] as 




(1.2) 



(cos 9) V'(O) + (sin 9) ip'iO) = 0, 



(1.3) 



^iV(0) + SiV'(0) = 0, 



(1.4) 
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such that the constant nx n matrices and Bi satisfy 

A^bI = B^AI (1.5) 

rank[^i Bi] = n, (1.6) 
i.e. AiBI is selfadjoint and the n x 2n matrix [Ai Bi] has rank n. 

Another formulation of the most general selfadjoint boundary condition at a; = for 
(1.1) was stated in [18-20] in terms of a constant n x n unitary matrix U2 as 

-By{0) + Al^p'{0)^0, (1.7) 

where the auxiliary constant nx n matrices A2 and B2 are given by 

A2:=^(C/2 + /n), B2:^'-{U2-In). (1-8) 

with In denoting the n x n identity matrix. It can directly be verified that the matrices 
A2 and B2 satisfy 

4^2=^24, BIB2=B2BI A2BI^B2AI 4^2 = 5+^2, (1.9) 
A2AI + B2BI = A2 + 1B2 = /n, A2 - 1B2 = U2. (1.10) 

We ourselves find it convenient to state the most general selfadjoint boundary condi- 
tion at X = for (1.1) in terms of constant n x n matrices A^ and B^ such that 

-sJV'(o) + 4V''(o) = o, (1.11) 

-5^^3 + 4^3 = 0, (1.12) 
4^3 + 53^3 > 0, (1.13) 

i.e. ^3-83 is selfadjoint and the selfadjoint matrix (^3^3 + B\B^) is positive. Note that 
(1.13) implies the existence of a unique positive matrix defined as 

E3:=(4^3 + 5JS3)^/^ (1.14) 
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such that is selfadjoint and invertible, and hence 

Es = El {EI)-\AIA3 + BIB3)E^' = In. 



(1.15) 



Let us define the matrices C3 and as follows 



B3 As 
A3 -B3 



H3 = C3 



E3' 



With the help of (1.14) and (1.15), it can be checked that H^H^ 
unitary. Thus, we must have H^H^ = l2n, yielding 



l2n and hence is 



A^Er^Al + B^E-^bI = /2„, B^E-'^Al - A^E-^B^ = 0. 



(1.16) 



We will mainly be working with the formulation given in (1.11)-(1.13) and hence later 
we will drop the subscripts in A3, B3, E3 and simply write A, B, E if there is no confusion. 

Let us note that one can multiply the boundary conditions stated in (1.4), (1.7), and 
(1.11) on the left by an invertible matrix D without changing the most general selfadjoint 
boundary condition at a; = 0. For example, for (1.11), by dropping the subscript 3, that 
left multiplication can be described via the transformation 

{A, B) ^ {A, B) {AD\ BD^), (1.17) 

for which we have 

{-B^A + A'^B) ^ -D-\-B^A + A^B){D^)-^, 

{A'^A + B^B) ^ L»-^(i^i + B^B){D^)-\ 

and hence (1.11)-(1.13) hold with {A,B) appearing instead of {A,B). Thus, the trans- 
formed pair {A, B) can be used instead of {A, B) in the formulation of the selfadjoint 
boundary condition. 

Our primary goal is to establish, under the most general selfadjoint boundary condition 
at a; = 0, the small- A; asymptotics of various quantities related to (1.1) such as scattering 
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solutions, the Jost matrix, the inverse of the Jost matrix, and the scattering matrix. The 
small-fc analysis for (1.1) has been lacking in the literature even though the relevant results 
are crucial in the study of the corresponding direct and inverse scattering problems. The 
direct scattering problem for (1.1) is to determine the scattering matrix and the bound- 
state information when the matrix potential V and the selfadjoint boundary condition is 
known. On the other hand, the inverse scattering problem is to recover the potential and 
the boundary condition from an appropriate set of scattering data. In some sense, our 
paper can be considered as a complement to the study by Agranovich and Marchenko 
[1], where the inverse scattering problem is analyzed only under the Dirichlet boundary 
condition but with attention to the behavior at /c = 0. Our study can also be considered as 
a complement to the study by Harmer [18-20] where the most general selfadjoint boundary 
condition (1.7) is used to investigate the inverse problem for (1.1) but the small-fc analysis 
is omitted. We refer the reader to [2,22] for similar small-/c analyses for the scalar radial 
Schrodinger equation, to [3,21] for the scalar full-line Schrodinger equation, to [4] for the 
matrix full-line Schrodinger equation, and to [5] for the radial Schrodinger equation with 
the most general selfadjoint boundary condition at the origin. 

Let us look at the definition of the Jost matrix J{k) given in (4.3). When V is 
selfadjoint and belongs to L^(R"'"), it is already known [1,4] that the right hand side in 
(4.3) is continuous in A; e C+, where we use C for the complex plane, 0+ for the upper half 
complex plane, and C+ := U R. Thus, J(0) exists. Let us also look at the definition 
of the scattering matrix S{k) given in (4.6) in terms of the Jost matrix J{k). In case J(0) 
is invertible, it is clear from (4.6) that S{0) = —In- However, if J(0) is not invertible, it 
is unclear whether S{k) is continuous at k = and what the value of S{0) is in case the 
continuity at A; = is assured. Our paper mainly concentrates on the case when J(0)~^ 
does not exist. We prove that S{k) is indeed continuous at A; = and we determine the 
value of ^"(0), which in general is different from —In- In case J(0) is invertible, our results 
reduce to the easy case with S{0) — —In- 
Let us note that J(0) is not invertible if and only if the determinant det[J(0)] is zero. 
In the scalar case (i.e. when n = 1) this is the analog of F0{O) = 0, where F^(A;) is the 
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Jost function appearing in (4.1). The case ^51(0) = is known as the "exceptional case," 
and the case F0{O) 7^ is known as the "generic case." Hence, in our paper we concentrate 
on the "exceptional case" for (1.1), namely the case when J(0) is not invertible. In the 
"generic case" it is already known and easy to see that S{k) is continuous at /c = and 
S{0) = —In- In the exceptional case, by expressing J{k) as in (6.17) in terms of a related 
matrix Z{k), and by writing the scattering matrix S{k) as in (6.19) in terms of Z{—k) and 
Z{k)~^, we are able to prove the continuity of S{k) at A; = and evaluate S{0). 

We remind the reader that the continuity of the scattering matrix in the exceptional 
case is not an easy matter. For example, in the full-line scalar case, Deift and Trubowitz 
[9] stated that the characterization of the scattering data given by Faddeev [12] might 
not hold and in fact even the continuity of the scattering matrix was not clear when the 
real- valued potential belonged to iv]^(R) and they introduced the stronger condition that 
the potential belonged to L2 (R) . The proof of the continuity of the scattering matrix when 
the potential belongs to I'i(R) was given later. For further details we refer the reader to 
[3,21] and the references therein. 

The matrix Schrodinger equation (1.1) has direct relevance to scattering in quan- 
tum mechanics involving particles of internal structures as spins, scattering on graphs 
[6,7,11,14,15,17,25-28], and quantum wires [23,24]. For example, the problem under study 
describes n connected very thin quantum wires forming a one-vertex graph with open 
ends. A linear boundary condition is imposed at the vertex and the behavior on each wire 
is governed by the Schrodinger operator. The problem has physical relevance to designing 
elementary gates in quantum computing and nanotubes for microscopic electronic devices, 
where, for example, strings of atoms may form a star-shaped graph. For the details we 
refer the reader to [23,24] and the references therein. 

Our paper is organized as follows. In Section 2 we show that the three selfadjoint 
boundary condition formulations given in Section 1 are equivalent. In Section 3 we in- 
troduce various n x n matrix solutions to (1.1) and state their properties relevant to the 
small-fc analysis of (1.1). In Section 4 we introduce the Jost matrix J{k) and the scatter- 
ing matrix S{k). In Section 5 we obtain various results that are crucial in determining the 



6 



small-A; asymptotics of the Jost matrix, its inverse, and the scattering matrix. In Section 6 
we provide the small-fc asymptotics for J(/c), J{k)~^, and S{k), and we prove that S{k) 
is continuous at k ~ 0. Finally, in Section 7 we provide some examples to illustrate the 
theory presented. 



2. EQUIVALENCE OF BOUNDARY CONDITION FORMULATIONS 

In Section 1 we have stated the three formulations of the most general selfadjoint 
boundary conditions at a; = for (1.1): 

(a) The formulation (1.4)-(1.6) stated in [23,24]. 

(b) The formulation (1.7) and (1.8) stated in [18-20]. 

(c) Our own formulation stated as (1.11)-(1.13). 

In this section we show that those three formulations are equivalent. 

Theorem 2.1 The three formulations (a), (h), and (c) of the most general selfadjoint 
boundary conditions at x — for (1-1) are all equivalent. 

PROOF: With the help of (1.17), we can relate (1.11)-(1.14) to (1.7) and (1.8) by letting 

U2 = {A^-iB^)E^\Al-iBl), 

and we can verify (1.7) and (1.8) with the help of (1.15) and (1.16). Hence, we have shown 
that (c) implies (b). Next, we will show that (b) implies (a). Let = —B\ and Bi = A\. 
Then (1.7) implies (1.4). Furthermore, the last equality in (1.9) yields (1.5). Let 



Co 



B2 A2 
A2 -B2 



We then get 



ClC2 



BIB2 + AIA2 BIA2-A\B2 

A\B2-BIA2 A\A2+BIB2 





'In 


■ 










where we have used (1.9) and (1.10). Thus, C2 is unitary and has rank 2n. As a result 
the block matrix [ s| AW has rank n. That matrix is in fact equal to [ —Ai B\ ] , and 



changing the signs in the first n columns does not affect its rank. Thus, [^i Bi] has 
rank n and (1.6) is satisfied. Hence, we have shown that (b) imphes (a). Finally, let us 
show that (a) implies (c). Let Bs — —A\ and A3 — Then, (1.4) implies (1-11), and 
(1.5) yields (1.12). Note that 



aIA3 + bIBs^[Bi Ai\ 



(2.1) 



and we need to show that the matrix product on the right in (2.1) is positive. This is 
indeed the case because that matrix product is itself a selfadjoint matrix and zero cannot 
be one of its eigenvalues. Otherwise, we would have a nonzero eigenvector v with the zero 
eigenvalue, implying 



= {v,[B^ A^ 





v} = { 


'Bf 




'bV 


A. 


A. 




A. 





(2.2) 



with (-, •) denoting the standard scalar product. However, (2.2) would then imply that 



Bl 



and hence the kernel of the matrix 



V. Consequently, the nullity of 
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would contain the nonzero vector 



B[ 



would be at least 1. Since the nullity and the rank 
must add up to n, the rank would have to be strictly less than n, violating the fact that 
the rank of that matrix is exactly n because of (1.6). Thus, (a) implies (c). | 

In the following proposition we state a fourth equivalent formulation of the most 
general selfadjoint boundary condition at a; = for (1.1). 

Proposition 2.2 The three formulations (a), (h), and (c) of the most general selfadjoint 
boundary condition are also equivalent to the formulation in terms of two constant n x n 
matrices A4 and S4 as 

-SlV(0) + 4V''(0) = 0, (2.3) 

such that the matrix C4 is unitary, where we have defined 

'B4 A4 ■ 
A4 -B4 



C4 



(2.4) 



PROOF: Because of Theorem 2.1, it is sufficient to prove the equivalence of (2.3)- (2.4) 
with (1.11)-(1.13). Suppose that (2.3)-(2.4) hold. Then, letting A3 = A4 and ^3 = B4, we 



get (1.11) and we obtain (1.12) and (1.13) from the unitarity of C4. Conversely, suppose 
(1.11)-(1.13) hold. By letting A4 = A^E^^ and S4 = B^E^^, we observe that (1.12) and 
(1.15) yield C4C4 = l2n, implying the unitarity of C4. | 

We have seen in (1.7) and (1.8) that the most general selfadjoint boundary condition 
for (1.1) can be stated in terms of a unitary matrix. There are certainly other choices for 
such a unitary matrix besides U2 appearing in (1.8). For example, in terms of a unitary 
matrix C/5, instead of (1.7) and (1.8) we can use 

-51^^(0) + 4^' (0) = 0, (2.5) 

where the auxiliary constant n x n matrices and B5 are given by 

A,:^'-(U,-Ul), B,:^^(U, + UI). (2.6) 

As seen from (2.6) we can simultaneously diagonalize into the form 

C/5 = diag{e^^i,e^^%...,e^^"}, (2.7) 

for some real valued parameters 9j. Then, the boundary condition given in (2.5) is separated 
into the n conditions given by 

(cos Oj) i;j{0) + (sin^,) V^^O) = 0, j = 1, . . . , n, 

where ipj denotes the jth column of the n x n matrix solution ip. Similarly, for the choice 
(1.8) for {A2, B2) in terms of a unitary matrix U2, hy diagonalizing U2 as in (2.7), we can 
express (1.7) as n separate boundary conditions given by 

[sin(^,/2)] V', (0) + [cos(^,/2)] V^^O) = 0, j = 1, . . . , n. 

3. PRELIMINARIES 

In this section we introduce certain n x n matrix solutions to (1.1) and state their 
properties that will be useful later on. We state the results without proofs and refer the 
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reader to the appropriate references such as [1,4] for details. Let us recall that we use the 
boundary conditions stated in (1.11)-(1.13) without the subscript 3. When V is selfadjoint 
and belongs to L{(R+), the matrix Schrodinger equation (1.1) has various n x n matrix 
solutions satisfying certain initial conditions or certain asymptotic conditions, and the 
existence of such solutions are already known. 



The Jost solution to (1.1) is the n x n matrix solution satisfying, for k e C+ \ {0}, 
the asymptotics 

f{k, x) = e^'^^iln + 0(1)], f'{k, x) = ik e^'^^iln + 0(1)], X ^ +00. (3.1) 
It satisfies the integral relation 

1 

f{k, x) = e'^^'In + Y dy sin k{y - x) V{y) f{k, y), 



X 



k 

and it is known [1,4] that f{k,x) and f'{k,x) are analytic in e C""" and continuous in 



k E C+ for each fixed x. The zero-energy Jost solution /(0,a:;) satisfies 

/"OO 

/(0,x) = /„+ / dy{y-x)V{y)f{0,y), 

J X 

and it is known [1,4] that /(0,a;) is a bounded solution to the n x n matrix-valued zero- 
energy Schrodinger equation 

-ip" + V{x)i; ^0, a;e(0,+oo), (3.2) 

satisfying 

/(O,x) = /n + o(l), /'(0,x) =o(l/x), X^+OO. (3.3) 
It is also known [1,4] that (3.2) has an n x n matrix solution g{0,x) satisfying 

g{0,x) = x[ln + o{l)], g'{0,x) = In + o{l), x^+oo. (3.4) 

Thus, the 2n columns of f{0,x) and g{0,x) form a fundamental set of solutions to (3.2), 
and any vector solution 0(a;) to (3.2) can be expressed as 

(l>{x) = f{0,x)^ + g{0,x)ri, x e {0,+oo), (3.5) 
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where the constant vectors ^ and rj in C" are uniquely determined by (f){x). We see from 
(3.4) and (3.5) that any solution to (3.2) that behaves as o{x) as a; — >■ +oo must be a 
bounded solution. 

There are various nxn matrix solutions to (1.1) defined via specifying some constant 
initial conditions at a finite x-value. As a result such solutions are analytic in k in the 
entire complex plane for each fixed x. Because of their analyticity such solutions are usually 
called "regular" solutions. The nxn regular solution (p{k, x) satisfies the initial conditions 

(p{k,0) = A, <p'{k,0) = B, (3.6) 

where A and B are the matrices appearing in (1.11). It satisfies the integral relation 

Sin JvCC X P"^ 

(p{k,x) = A cos kx + B — h — / dy sin k{x — y)V{y) (fi{k,y). (3-7) 

k k Jq 

Let us define two additional regular nxn matrix- valued solutions to (1.1), namely 
C{k, x) and S{k, x) with initial conditions at a; = o, at which the matrix /(O, o) is invertible. 
The existence of such an a-value is assured by the fact that /(O, x) = + o(l) as x — ?■ +oo 
and hence /(0,a;) is invertible at least for large x-values. In fact, if /(0,a)~^ exists, then 
we must have the existence of f{k, a)~^ in the vicinity of A; = in C+. This is because for 
each fixed x-value it is known [4] that f{k,x) is a continuous function of A; e C+. Hence, 
det[/(/c, a)] is a continuous function of k and if it is nonzero at A; = it must be nonzero 
in the vicinity of k — 0. Thus, we conclude that 

f{k, a) = /(O, a) + o(l), /(/c, a)-i = /(O, a)-^ + o(l), /c ^ in C+. (3.8) 

The cosine-like solution C{k,x) satisfies the initial conditions 

C{k, a) = In, C'{k, a) = 0, (3.9) 
and the sine-like solution S{k,x) satisfies 



S(k,a) = 0, S'{k,a) = In. 
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(3.10) 



Thus, we have the integral representations 

1 r 

C{k,x) = I cos k{x - a) + - dy sin k(x - y)V(y) C{k,y), (3.11) 

^ J a 

S{k, x)=I +1 T'^y - ^) ^^^) ^3.12) 

Note that we suppress the dependence on a in our notation for such solutions. 

We define another n x n regular solution to (1.1), uj{k,x), which satisfies the initial 
conditions 

uj{k, a) = /(O, a), uj'{k, a) = f'{0, a). (3.13) 
Again we suppress the dependence on a in our notation for oj{k, x). Note that 

a;(0, x) = /(O, x), xe R+, (3.14) 

because both sides satisfy (1.1) when k = and they both satisfy the same initial conditions 
at X = a given in (3.13). It is seen from (3.9), (3.10), and (3.13) that 

u{k, x) = C{k, x) /(O, a) + S{k, x) /'(O, a), (3.15) 

where f{k,x) is the Jost solution appearing in (3.1). 

Let us note that our regular solutions satisfy for /c e C 

(f{—k,x) = (fi{k,x), C{—k,x) = C{k,x), S{—k,x) = S{k,x), u){—k,x) = oj{k,x). 

(3.16) 

This is because k appears as k'^ in (1.1) and the initial values of those solutions are 
independent of k, as seen from (3.6), (3.9), (3.10), and (3.13). 

Associated with (1.1) we have the adjoint equation 

-iP^" + il;^V{x) = {k*fiP'^, xe{0,+oo), (3.17) 

where we have used (1.2) and an asterisk denotes complex conjugation. Note that if '0(/c, x) 
is any solution to (1.1), then ip{±k*,x)'^ is a solution to (3.17). Let us also add that if 
'ip{k,x) has an analytic extension from k e H to k e C^, then 'ip{—k,x)^ has also an 
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analytic extension from A; e R to A; e C""", and in fact that extension becomes equal to 
ijj{—k*,x)^ for k e C^. A consequence of this is the following. Since it is already known 

that f{k,x) and f'{k,x) are analytic in A; G C+, /(— and /'(— have analytic 
extensions from /c e R to /c e C"*" given by f{—k*,xy and f {—k* , , respectively. 

Let [F; G] := FG' — F'G denote the Wronskian. It can directly be verified that 
for any n x p solution ip{k,x) and any n x q solution (/){k,x) to (1.1), the Wronskians 
[(j){k* , x)"^ ; tp{kj x)] and [(/){— k* ^ x)"^ ; tp{k, x)] are both independent of x. By evaluating the 
values of the Wronskians at a; = and x = +00, we can obtain various useful identities. 
For example, we have 

[f{±k, x)^; f{±k, x)] = ±2ikln, /c e R, (3.18) 
[/(-r,x)+;/(A;,x)] = 0, keC+. (3.19) 



4. THE JOST MATRIX AND THE SCATTERING MATRIX 

In this section we introduce the Jost matrix and the scattering matrix for (1.1) with a 
selfadjoint matrix potential V in L]^(R+) with the selfadjoint boundary condition (1.11)- 
(1.13). We also present certain preliminary results needed later on to analyze the small-A; 
limits of these two matrices and of the inverse of the Jost matrix. 

Recall that the Jost function Fe corresponding to (1.3) in the scalar case, i.e. when 
n = 1 in (1.1), is defined with the help of the Jost solution f{k,x) as [5,16,29,31] 

r -i[/'(/c,0) + (cot^)/(A;,0)], ^e(0,7r), 
Fe{k) := <^ (4.1) 

We will define the matrix analog of the Jost function, which is called the Jost matrix, so 
that it reduces to the familiar Jost function when n — 1. Recall also that the scattering 
matrix in the scalar case is defined as [5,16,29,31] 

(4.2) 

F0{-k) 



Se{k) := 



I Feik) ' 
13 



6» = 7r. 



The reason behind the sign difference in (4.2) in the Dirichlet case (i.e. when ^ = tt) is 
that (4.2) ensures that 5'(/c)— >lasF— >0, which is a consequence of the fact that the 
perturbed and unperturbed Hamiltonians satisfy the same selfadjoint boundary condition 
at X = 0. We will define the scattering matrix by generalizing (4.2) to the matrix case. For 
simplicity, we will suppress the dependence of the Jost matrix and the scattering matrix on 
the boundary-condition parametrization (A, B), and we will use the notation J{k) for the 
Jost matrix instead of J{A,B){k) and also write S{k) for the scattering matrix instead of 
S(^A,B){k). Note that we earlier used S{k, x) in (3.9) to denote the sine-like regular solution 
to (1.1), which should not be confused with the notation S{k) used for the scattering 
matrix. 

Define the Jost matrix J{k) for k E C+ as 

J{k) := [/(-r , x)t; ^{k, x)] = /(-r , 0)ti? - f'{-k, 0)U, (4.3) 

where /(/c, x) is the Jost solution appearing in (3.1), (fi{k, x) is the regular solution appear- 
ing in (3.6), and A and B are the matrices appearing in (1.11)-(1.13) and (3.6). Note that 
J is not uniquely determined by the potential V and the selfadjoint boundary condition 
(1.11)-(1.13). This is because (1.11)-(1.13) are invariant under the transformation (1.17), 
and hence we have J i— )■ JD"^ under (1-17) indicating that the definition for J in (4.3) is 
unique up to a right multiplication by a constant invertible matrix. On the other hand, 
such a postmultiplication does not change the zeros in C"*" of the determinant of J{k). 
Those zeros correspond [1,19] to the bound-state energies of (1.1) with the boundary con- 
dition (1.11)-(1.13), and hence the bound state energies are still uniquely determined by 
(4.3). 

Theorem 4.1 If V is selfadjoint and belongs to L\(R,'^), then the Jost matrix J{k) is 
invertible for A; e R \ {0}. 

PROOF: Even though a proof is available [19], for the benefit of the reader we outline a 
proof of our own. For A; e R define 



L{k) := f'{-k, 0)^BE-^ + f{-k, 0)Ue-2, 
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(4.4) 



where E is the matrix E3 appearing in (1.14). With the help of (1.16), (3.18), (4.3), and 
(4.4) one can show that 

J{k) L{k)^ - L{k) J{ky - [/(-fc, x)^- f{-k, x)] = -2ikln, A; e R. (4.5) 

If J{k) were noninvertible at some real nonzero fco, then the rows of J(/co) would be linearly- 
dependent and hence we would have tt^ J(/co) = for some nonzero vector u G C"^ as well 
as J(A;o)^w = 0. However, because of (4.5) this would imply 

= u^J{ko) Liko^u - u^L{ko) Jiko^u = -2ikou^u = -2iko \ \u\f ^ 0, 

which is a contradiction. Thus, J{ko) must be invertible. | 
The scattering matrix S{k) is defined as [18-20] 

S{k) := -Ji-k) J{k)-^, keR\ {0}, (4.6) 

and it is uniquely determined by the boundary condition and the potential V. Even though 
J{k) is uniquely defined only up to a right multiplication by a constant invertible matrix, 
the unique determination of S{k) is assured because S{k) remains invariant under the 
transformation (1.17). Note that the domain of J{k) is A; e C+ because /(— A;,0)^ and 
f'{—k, 0)''" have analytic extensions from k E H to k e 0+ and the values of those exten- 
sions are f{—k*,0)^ and f'{—k*,0)^, respectively. On the other hand, in general S{k) is 
defined only for real k because J{—k) in general cannot be extended from k eUtok e C^. 
Furthermore, the existence of S{k) when k = needs to be studied separately because, as 
we have seen in Theorem 4.1, the existence of J{k)~^ is assured only for k e R\ {0} and it 
cannot easily be inferred from (4.6) whether S{k) has a limit as A; — > when J(0)~-^ does 
not exist. 

In order to understand the small-/c behavior of J{k) and S{k), it is instructive to 
analyze first the case when the potential V is identically zero in (1.1). In that case, we 
have f{k,x) = e^^^In, and hence (4.3) and (4.6) yield 



J{k) = B-ikA, [J{k)]-^ = {B-ikA)-'^, Sik) = -{B + ikA)(B-ikA)-K 
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Let us use the representation (2.6) for {A,B) with the diagonal form of U given in (2.7). 
We then obtain 

A = — diag{sin^i, . . . ,sin^„}, B = diag{cos^i, . . . ,cos^n}5 
J{k) = diag{ Ji(/c), . . . , Jn{k)}, S{k) = diag{5i(A;), . . . , Sn{k)}, 
where we have defined 

As seen from (4.7), in the Dirichlet case (i.e. when 9j = n) we have 

Mk)^-i, [j,(fc)]-i = -i, s,{k)^-i. 

On the other hand, in the Neumann case (i.e. when dj = 7r/2) we have 

jj{k) = ik, [jjm-' = ^, sj{k) = i. 

Note that, in the Neumann case, Jj{0) vanishes hnearly as ^ and it is not an invertible 
matrix; however, Sj{0) is still well defined because J{—k)[J{k)]~^ has a well-defined limit 
as A; — > 0. It is somehow disturbing that in the Dirichlet case, Sj{0) ^ 1 and in fact 
'S'j(O) = ~1) which is exactly the opposite of the scalar case as seen from (4.2). The 
explanation for the discrepancy is that the unperturbed Hamiltonian in the matrix case 
is chosen to satisfy the Neumann boundary condition, which is compatible with the time- 
dependent derivation of the scattering matrix and motivated by applications in quantum 
wires; for further elaboration on this point we refer the reader to p. 1566 of [24]. 

5. SMALL /c-BEHAVIOR 

In preparation for the analysis of the small-energy behavior of the Jost matrix J{k), 
its inverse J{k)~^, and the scattering matrix S{k)^ in this section we establish the small- A; 
asymptotics of various quantities related to the regular solutions to (1.1). 

We are interested in analyzing the Jost matrix J{k) as /c ^ in C+. From (4.3) we 
see that 

J(0) = /(0,0)tS-/'(0,0)U, (5.1) 
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and we would like to determine how fast J{k) approaches J(0) and whether J{k) ^ exists 



at = and determine its behavior as — > from C+. We would like to know about such 
small-/c behaviors when V is selfadjoint and belongs to L\{'R,'^). 



As stated before (3.8), f(k, a) is invertible in the vicinity of A; = in C+ for some a 
value. In (4.3) we have defined the Jost matrix in terms of a Wronskian whose value is 
independent of x. As we see below we can write J(fc)^ in terms of Wronskians evaluated 
at X = a and involving the solutions f{k,x), (p{k,x), and oj{k,x) appearing in (3.1), (3.6), 
and (3.13), respectively. 

The following result will be needed later on. By a generic constant, we mean a constant 
that does not necessarily have the same value in difi'erent appearances. 

Proposition 5.1 If V is selfadjoint and belongs to L];(R"''), then the regular solution 
u){k,x) to (1-1) appearing in (3.13) satisfies 



\uj{k,x)-Uj{<d,x)\\<c( *^|[^^ ) ^{lra[k]){x-a) ^ k G C+ , X > U, (5.2 

y 1 -|- I A; I {^x cf') / 



where c is a generic constant. 

PROOF: From (3.11), (3.12), and (3.15) we have 

uj{k, x) = /(O, a) cos k{x - a) + /'(O, a) ^ ^ / ~ ^^^^ ^^^^ 

(5.3) 

Note that (5.3) yields 

a;(0, x) = /(O, a) + {x- a) f'{0, a) + [ dy {x - y) V{y) uj{0, y), (5.4) 

J a 

and u){0,x) = f{0,x) by (3.14). Thus, from (5.4) and its x-derivative, with the help of 
f'{0,x) = o{l/x) as a; — >■ +oo, we obtain 

/•OO 

/ dyV{y)uj{0,y) = -f{0,a), (5.5) 

J a 

r-OO 

/ dyy V{y) oo{0, y) ^f{0,a)-a /'(O, a) - (5.6) 

J a 

17 



Let us write (5.5) as 

j-X nOO 

/'(0,a) = -/ dyV{y)u{0,y)- dyV{y) u{0,y). 

J a J X 

Using (5.3), (5.4), and (5.7), we get 

uj{k, x) - uj{0, x) = Ki + K2 + Ks + K4, 

where we have defined 



(5.7) 



Ko := 



Ki := /(O, a) [cos k{x — a) — 1] , 

sin k{x — a) 
k{x — a) 



(5.8) 

(5.9) 
(5.10) 



(a; -a) / dyV {y) uj{Q,y), 

J X 

1 r 

:= — / dy [sinkix — y) — kix — y) — sinkix — a) + kix — a)\V{y)u{Q,y), (5.11) 

k J a 

K^-=\j^ dy[suik{x-y)]V{y)[u{k,y)-u{Q,y)]. (5.12) 



For 2 e C+ we have 
I sin^l < 



c 


z 


glm[z] 


1 + 


Z 





sm^ 



< 



c 


z 


2 glm[z] 


(1+ 


z\Y 



11 — COS < 



c 


z 


2 glm[z] 


(1 + 


A? 



(5.13) 



Using (5.9), the third estimate of (5.13), and the fact that f{0,a) is bounded, we get 



{l + \k\{x-a)f 



k e C+, x >a. 



(5.14) 



Note that 



(a; -a) / dyV{y)u{0,y) 



< 



< 



dy {y - a) V{y)uj{0,y) 
dyyV{y)uj{0,y) 



(5.15) 



The norms in (5.15) are bounded by a constant due to the facts that V e Lj(R+) and 

a;(0, x) is bounded as a result of (3.3) and (3.14). Thus, from (5.10), (5.15), and the second 
estimate in (5.13), we get 

c|fc|2(a;-a)2e(i"^W)(^-'*) 



{1 + \k\{x - a))^ 
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k e C+, X > a. 



(5.16) 



Let us now estimate when k e C+ and x > a. We write (5.11) as 

nX nX — a 

K3= dy dz [1 - coskz] V{y)u{0,y), 

J a J x—y 



V 

2 //'I I „\2 



and use the third estimate of (5.13) and the fact that x i-> x /{I + x) is an increasing 
function of x when a; > 0, to obtain 

11^311 < '(U|4-ay L <'"lin.)IIIMO,.)l|. (5.17) 
Since w(0,y) is bounded and V e Lj;(R+), from (5.17) we obtain 



for a generic constant c. Let us now estimate when A; e C+ and x > a. Letting 

C{k,x) := e-(^™['=l)^^-'')||a;(A;,a;) -a;(0,a;)||, (5.19) 

from (5.12) we get 

g-(im[fe])(a.-a)||^^|| < 1 T g_(im[fc])(a.-y) [gin A;(a; - y) | | |F(y) 1 1 C(/c, 2/). (5.20) 

FI Ja 

Using the first estimate of (5.13) in (5.20) and the fact that x i->- x/{l + x) is an increasing 
function of x when a; > 0, we obtain 

1 I < iyymy)\m,y). (5.21) 

Using (5.14), (5.16), (5.18), and (5.21) in (5.8) we get 

X c\kP(x — a)^ c \k\(x — a) , ... . ,^ 



Setting 



, [1 + \k\ix — a)Y ^ , , , 

^ - |J|2 & ((^^ ^ ' 5.23 

c \k\^{x — ay 



we can write (5.22) as 

IS/,, „\2 
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From (5.24) we obtain 

Xik,x)<l + c[ dyy\\Viy)\\xik,y), (5.25) 

where we have used, for < a < y < x, the estimate 

|/cp(y — a)^ ^ \k\{x ~ a) \k\{y — a) ^ \k\{x — a) , , 

(1 + \k\{y - a))2 - (1 + \k\{x - a)) (1 + \k\{y - a)) " (1 + \k\{x - a)) ' ' ^' 

based on the fact that x ^ x/{l + x) is an increasing function of x when x > 0. Applying 
GronwaU's lemma to (5.25) and using the fact that V G L\{K+), we obtain x{k, x) < c for 
some generic constant c, which is not necessarily equal to the generic constant c in (5.25). 
Thus, using (5.19) and (5.23) in x{k,x) < c we obtain (5.2). | 

Let us define 

P{k):=[u{-k*,x)^-J{k,x)i (5.26) 

where we note that the Wronskian in (5.26) is independent of x, and hence with the help 
of (3.13) and (3.16) by evaluating that Wronskian at a; = a we get 

Pik) = /(O, a)t/'(A;, a) - /'(O, a)t/(A;, a). (5.27) 

Note that P{k) has an analytic extension from k e H to k e C"*" because f{k, a) and 
f'{k, a) possess that property as well. It is difficult to obtain useful information from 
(5.27) as /c ^ because for V G Li(R+) we can only say that 



f{k, x) = /(O, x) + o(l), f{k, x) = /(O, x) + o(l), /c ^ in C+. (5.28) 

In the proposition below we evaluate the small-/c asympotics of P{k) by evaluating the 
Wronskian in (5.26) at a; = +oo. This result will be useful in evaluating the small-/c limit 
of the Jost matrix J{k). 

Proposition 5.2 If V is selfadjoint and belongs to L}(R+), then the matrix P{k) given 
in (5.26) satisfies 



P{k) = ikin + o{k), /c ^ in C+. (5.29) 
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PROOF: We can evaluate the asymptotics of uj{—k*,x) and u'{—k*,x) as x ^ +00 from 
(5.3). Furthermore, the asymptotics of f{k,x) and f'{k,x) as x ^ +00 are available from 
(3.1). Using those asymptotics in (5.26) we obtain 

/•CO 

P(/c) = iA;e^'="/(0,a)t-e^^"/'(0,a)+- / dy e^''^ uj{-k* ,y)^ V{y). (5.30) 

J a 

Let us break the right hand side in (5.30) into three terms and write 

Pik) = Pi{k) + P2{k) + Ps{k), (5.31) 

where we have defined 

Piik) := ike''"' /(O, a)+ - e''''' /'(O, a)\ (5.32) 

poo 

P2{k):=- / dye"'y[uj{-k*,y)^ -uj{0,yy]V{y), (5.33) 

J a 

P3{k):=- / dye"'yuj{0,y)^V{y). (5.34) 

J a 

From (5.2) and the Lebesgue dominated convergence theorem, it follows that P2{k) = o{k) 



as A; ^ in C+ whenever V e Lj(R+). Using (5.5) and (5.6) in (5.34) we have 

Pm = /'(O, a)t + ik[In + a /'(O, a)t - /(O, a)t] + P^k), (5.35) 
where we have defined 

POO 

P^{k) ■■=- dy [e'''y - 1 - iky] a;(0, y)^ V{y). (5.36) 

J a 



Note that for any 2; e C+ we have 

I 1 2 

\e'' - 1\ < c\z\, le'^ -l-iz\< (5.37) 
where c denotes a generic nonnegative constant independent of the complex number z, and 



hence the second inequality in (5.37) helps us to get P4{k) = o{k) as A; —> in C+. Thus, 
using (5.31)-(5.36), we have 

P{k) = ikin + [ike''''' - ik] /(O, a)^ + [In + ika - e'''"] /'(O, a)^ + o(/c), 
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and because of (5.37) each of the coefficients of f{0,ay and f'{0,a)^ is 0{k'^) as A; in 
C+. I 

Even though (5.28) does not provide much extraordinary information, the following 
theorem shows that f'{k^x)[f{k,x)]~^ is differentiable at A; = at any fixed x value 
where the matrix /(0,a;) is invertible. The results stated in the next theorem are the 
generalization to the matrix case of similar results in the scalar case [2,3]. 

Theorem 5.3 Assume that V in (1.1) is self adjoint and belongs to L];(R+). If the constant 
matrix /(0,a) is invertible, where f{k,x) is the Jost solution appearing in (3.1), then, 
f'{k,a) [f{k,a)]~'^ is differentiable at k = and we have 

f'{k, a) f{k, a)-' = /'(O, a) /(O, a)-' + ik[f{0, a)" ^/(O, a)-' + o{k), A; ^ in C+. 

(5.38) 

If instead, f'{0, a) is invertible at some a e R""", then the Jost solutions satisfy 

f{k, a) f'{k, a)-' = /(O, a) /'(O, a)-' - ik[f'{0, a)-']^f'{0, a)-' + o{k), A: ^ in C+. 

(5.39) 

PROOF: In case f{0,a) is invertible, by the continuity of the determinant of f{k,a), we 
must have f{k, a) invertible in C+ in the vicinity of A; = 0. Thus, from (5.27) we get 

fik, a) fik, a)-' = [/(O, a)t]- V'(0, a)t + [/(O, a)^]-'P{k) f{k, a)-\ (5.40) 

Note that (3.19) implies that 

[/(O, a)t]- V'(0, a)t = /'(O, a) /(O, a)-\ (5.41) 

Applying (5.29) in (5.40) and using (3.8) and (5.41), we obtain (5.38). In a similar way, 
(5.39) is proved. | 

Next, we express the Jost matrix J{k) defined in (4.3) in terms of the Jost solution 
f{k,x), the regular solutions Lp{k,x), the regular solution u;{k,x), and the matrix P{k) 
appearing in (3.1), (3.6), (3.13), and (5.26), respectively. 

Proposition 5.4 Assume that V in (1.1) is selfadjoint and belongs to L\(R'^). The Jost 
matrix can be written as 

J{k) = Ti{k) + T2{k), ke C+, (5.42) 
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where we have defined 

T,{k) := -P(-r)V(0,a)-V(fc,a), (5-43) 
T^ik) := /(-r , a)t[/(0, a)-i]t[a;(-r , x)t; (^(/c, a;)], (5.44) 

and recall that the value of the Wronskian appearing in (5.44) independent of x. 
PROOF: Using (4.3) we can write J{k) in terms of Wronskians evaluated at a; = a to have 

Jik) =[fi-k\x)^;^ik,x)]l^^ 

=/(-r , a) V(fc, a) - /'(-r , a) a) 
=T^{k) + T2{k)+T3{k), 

where we have defined 

Ti(fc) - (/'(-r , a) V(0, a) - f{-k*, a) V'(0, a)) /(O, a)" V(^, «), (5-45) 

T2(fc) := /(-r , a)t[/(0, a)-i]t (/(O, a) V(fc, a) - /'(O, a)tvp(fc, a)) , (5.46) 

nik) := /(-r , a)t[/(0, a)-i]t (/'(O, a) V(0, a) - /(O, a)t/'(0, a)) /(O, a)" V(A;, a). 

(5.47) 

Using (5.27) in (5.45) we see that Ti(k) in (5.45) can equivalently be written as in (5.43). 
Using (3.13) in (5.46), we see that T2{k) in (5.46) can equivalently be written as 

T2{k) =f{-k\ a)t[/(0, a)-i]^ {uj{-k\ a) V(fc, a) - u'{-k*, a) V(A:, a)) 
=/(-r , a)t [/(O, a)- V[^(-A:*, a;)t; ^(/c, x)] |^^^, 

which is equivalent to (5.44) because the Wronskian in (5.44) is independent of x and can 
be evaluated at x = a. Finally, using (3.19) with k — in (5.47) we see that T2{k) = 0. | 

In the next theorem we evaluate the Wronskian [uj{—k*,xy;(p{k,x)] appearing in 
(5.44). 

Proposition 5.5 Assume that V in (1.1) is self adjoint and belongs to L];(R+). Then, the 
Wronskian appearing in (5.44) has the small-k asymptotics 

[a;(-r,a;)^;(^(/c,a;)] = J(0) + O(/c^), /c ^ in C+, (5.48) 
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where J{k) is the Jost matrix defined in (4-3). 

PROOF: Since the value of the Wronskian in (5.48) is independent of x, we wiU evaluate 
its value at a; = 0. By writing 

*,x)^ - oo{0, x)^; ip{k, x)] + [a;(0, a;)+; ^(/c, x% (5.49) 

from (3.6), (3.14), and (5.1) we see that the second Wronskian on the right side in (5.49), 
when a; = 0, yields 

[a;(0,x)t;^(A;,x)]|^^^ = J(0). 

Next, we evaluate at x = the value of the first Wronskian on the right side in (5.49). 
The x-derivative of that first Wronskian, with the help of (1.1) and (3.17), can be directly 
evaluated as 

^[oj{-k\x)^ -u:{^,x)^;ip{k,x)] = k'^ u{^,x)^ ip{k,x). (5.50) 
Integrating (5.50) over the interval [0, a], and then using (3.13) and (3.14), we obtain 

[a;(-r , x)t - a;(0, x)^- ^{k, x)] = -k^ [ dy /(O, <f{k, y). 

Jo 

Thus, we have the estimate in (5.48). | 

Using (5.4), (5.29), (5.42)-(5.44), and (5.48) we have the following conclusion. 

Corollary 5.6 Assume that V in (1.1) is selfadjoint and belongs to L^(R+). Then, the 
Jost matrix J{k) appearing in (4-3) has the small-k behavior 

f{0, a)^[/(-/c*, a)^]-^ J(/c) = J(0) - ik /(O, a)"^ v?(0, a) + o{k), /c ^ in C+, (5.51) 

where f{k,x) and (f{k,x) are the Jost solution and the regular solution appearing in (3.1) 
and (3.6), respectively, and a is any point where the matrix /(O, a) is invertible. 

In order to study the small-A; limit of J{k)~^., we will next concentrate on the 0{k)- 
term appearing in (5.51), namely /(O, a)~^ <^(0, a). We use Ker [J(0)] to denote the kernel 
of the matrix J(0). 

Proposition 5.6 Assume that V in (1.1) is selfadjoint and belongs to L{(R+). Then, the 
following are equivalent: 
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(a) The vector u e is an eigenvector of the zero-energy Jost matrix J(0) with the zero 
eigenvalue, i.e. u & Ker[J{0)]. 

(b) (f'{0,-\-oo)u = 0, where (f{k,x) is the regular solution to (1.1) appearing in (3.6). 

(c) (/?(0, x) u is hounded for x G [0, +oo). 
PROOF: From (4.3) we see that 

J(0) = /(O, x)^ ^'(0, x) - /'(O, x)t ^(0, x), (5.52) 

where the quantity on the right side in (5.52) is independent of x. From (3.5) it follows that 
each column of </7(0,a;) is a linear combination of columns of /(0,a;) and g{Q,x). Hence, 
there exist constant n x n matrices a and P such that 

ip{0,x) = f{0,x)a + g{0,x) f3, a; e R+. (5.53) 

From the x-derivative of (5.53) we get 

ip'{0,x) ^ f'{0,x)a + g'{0,x) f3, x e R+. (5.54) 

Using (3.3) and (3.4) in (5.54) we see that 

(/?'(0, +oo) = /3. (5.55) 

Inserting (5.53) and (5.54) on the right side of (5.52), evaluating the resulting expression 
as a; ^ +oo, and using (3.3) and (3.4), we obtain 

J(0) = /3. (5.56) 

Thus, from (5.55) and (5.56) we conclude that 

J(0) = ^'(0,+oo), (5.57) 

and hence the equivalence of (a) and (b) are established. Note that, from (3.3), (3.4), 
and (5.53) it follows that (p{0,x)u is bounded if and only if ^tt = 0, which happens if 
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and only if (fi'{0,+oo)u = as a result of (5.55). Thus, the equivalence of (b) and (c) is 
established. | 

Let us note that we can express </7'(0, +00) in (5.57) in another form. Letting k ^ 
in (3.7) we get 

ip{0,x) = A + Bx+ f dy{x-y)V{y)if{Q,y), (5.58) 
and from the x-derivative of (5.58) we have 

ip'{0,x) = B+ [ dyV{y)ifi{0,y), (5.59) 

We know from (3.3)-(3.5) that (p{0,x) can grow at most as 0{x) as x ^ +00 and hence 
the integral in (5.59) exists as a; — >■ +00, and from (5.57) and (5.59) we get 



J{0) = cp'{0,+oo) = B+ / dyV{y)ip{0,y). 

Jo 



Proposition 5.7 Assume that V in (1.1) is selfadjoint and belongs to Lj;(R+). Then, for 
any vector u in Ker[J{0)] there exists a unique vector^ in Ker[J{0)'^] such that 

if{0,x)u^f{0,x)t (5.60) 

The map from Ker[J{0)] to Ker[J{Oy] is a bijection. 

PROOF: By expressing (fi{0,x)u as in (3.5) we get 

ip{0, x)u = /(O, x)^ + t/(0, x) r), 

and hence (f{0,x)u is bounded, equivalently as stated in Proposition 5.6, u e Ker [J(0)], 
if and only if = 0. Thus, the mapping tt i-> ^ is identified with (5.44). Let us now show 
that ^ e Ker [J(0)t]. Using (5.1), we have 

J(0)te = [5t/(0, 0) - A^f'iO, 0)]e. (5.61) 

On the other hand, from (5.44) and its derivative, with the help of (3.6), we see that 



/(0, 0) e = ^{0, 0)u = Au, /'(O, 0) e = 0)u = Bu. (5.62) 

26 



Using (5.62) in (5.61) and imposing (1.12) we get J{0)^^ = 0. Thus, ^ e Ker [J(0)+]. Note 
that the map it i->- ^ is a hnear map from Ker [J(0)] into Ker [J(0)^] because, as seen from 
(5.60), we have 

^ = /(0,a)"V(0,a)u, « e Ker [J(0)]. (5.63) 

The map tt i-)- ^ from Ker [J(0)] to Ker [J(O)'''] is one-to-one because it has zero kernel as 
seen by the foUowing argument. If ^ = in (5.60), then we must have f{0, x) u = and 
ip'{0,x)u = 0. In particular, at a; = with the help of (3.6) we then have Au = and 
Bu = 0. We in turn get A^Au = and B'^Bu = 0, and hence {A^A + B'^B)u = 0, which 
yields u = because of (1.13). Furthermore, Ker [J(0)] and Ker [J(O)'''] have the same 
dimension, and hence the map tt i-> ^ from Ker [J(0)] to Ker [J(0)^] is a bijection. | 

6. SMALL /c-BEHAVIOR OF J{k)-^ AND OF S{k) 

In this section we establish the small-/s asymptotics of the Jost matrix J{k), its inverse 
J{k)~^, and the scattering matrix S{k). As we will see, J{k) is continuous at = 0, J{k)~^ 
has an 0{l/k) singularity at = if J(0) has a zero eigenvalue, J{k) is continuous at 
/c = if zero is not an eigenvalue of J(0), and that S{k) is continuous at k = whether or 
not zero is an eigenvalue of J(0). 

In order to analyze the small-/c behavior of J{k)~^, we will analyze (5.51). Let us 
write (5.51) as 

F{k) = J{0) -ikR + o{k), A; ^ in C+, (6.1) 

where we have defined 

Fik) := /(O, a)t[/(-r , a)-i]t J(A;), R := /(O, a)'' <^(0, a). (6.2) 

We will equivalently analyze the behavior of F{k)~^ as A; — > e C+. As we have seen in 
Proposition 5.7, the restriction of R to Ker [J(0)] yields an invertible map. Among the n 
eigenvalues of J(0), let us assume that the zero eigenvalue has geometric multiplicity ^ with 
a possibly larger algebraic multiplicity v. In other words, J(0) has /i linearly independent 
eigenvectors corresponding to the zero eigenvalue. Thus, J(0) d = if and only if e C" 
is such an eigenvector, in which case we have Rv ^ 0. 
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Let us choose a Jordan basis [10] for the matrix J(0) as foUows. Assume that there 
are k Jordan chains and hence k blocks in the Jordan canonical form of J(0). Let us use 

the index a for a = 1, . . . , k to identify the Jordan chains and assume that the ath chain 
consists of Ua vectors Uf^j for j = 1, . . . , Let us use Aq, to denote the eigenvalue of J(0) 
associated with the ath chain, where the eigenvalues may be repeated. We have 

[J(0) - Xa] Ual = 0, 

[J(0) - Xa] Uaj = Ua(j-l), j = 2, . . . , TT-q, 

and hence Uai is an eigenvector and Uaj for j = 2, . . . , Ua are the generalized eigenvectors. 

Since we assume the zero eigenvalue has geometric multiplicity fi, without loss of 
generality we let Aq. = for a = 1, . . . , /x and Aq, 7^ for a = /x + 1, . . . , We order the 

vectors in the Jordan basis according to the rule that Uaj comes before tt^g if and only if 
ct < /3 or a = /3 and j < s. Thus, {un, U21, . . . , w^i} forms a basis for Ker [J(0)], and our 
Jordan basis is given by the ordered set {uaj}, i.e. 

{uii,Ui2, Wlni, "21, '"22, • • • , M2n2 , • • • '"«2, • • • , W«n,}- (6.4) 

The corresponding adjoint Jordan basis {vaj} satisfies vl,jUpt = SapSjt, with 5jt 
denoting the Kronecker delta, and the indices a and p referring to the Jordan blocks. The 
vectors Vaj satisfy 

[J(0)t-A*]^Q.„ =0, ^g^^ 

[J(0)+ -X%]Vaj ^Ua(j+i), j = l,...,nQ - 1. 

Thus, {vini 5 ■J^2n2 5 • • • j'^m^m) ^rms a basis for Ker [J(0)'l']. The adjoint Jordan basis is the 
ordered set {vaj}, i-e. 

it'll, t'12, • • ■,Vinj_,V2l,V22, ■ ■ • , t'2n2 , • • •^k1,^«2, • • • j^Kn^}- (6.6) 

Let S denote the matrix whose columns are given by the elements of the ordered set 
{uaj} in (6.4). Then, is exactly the matrix whose rows are given by the elements of 
the ordered set {f^j}, with the ordering given in (6.6). Thus, the Jordan canonical form 
of J(0) is given by 

5-^J(0)5 = e5=iJ„jAQ), (6.7) 
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where Jnc('^a) is the ria x Jordan block with Aq, appearing in the diagonal entries and 
one in the superdiagonal entries. Since the first /x Jordan blocks Jn„(-^a) are associated 
with the zero eigenvalue and the remaining {n — ^) blocks are associated with nonzero 
eigenvalues, each Jnc{}^a) is an na x Hq, matrix given by 

1 ... 0- 



1 















1 

An- 





1 


















1 
0. 







(6.8) 







A. 



CK = + 1, . . . , K. 



(6.9) 



Let us use a tilde to denote the transformation via iS, i.e. M := S ^MS for any nxn 
matrix M. Let us apply this transformation on the matrix F{k) appearing in (6.1) and 
(6.2). Then (6.1) yields 



F{k) = J{0) -ikR + o{k), k^OmC+. 



(6.10) 



By inspecting (6.7)- (6. 10) we see that there are exactly /j, columns of F{k) behaving as 
0{k) as k ^ and each of the remaining (n — /i) column vectors contains at least one 
entry that has a nonzero limit as — > 0. 

Our next goal is to move all the entries with 1 appearing in the superdiagonal in the 
first iJ, Jordan blocks in (6.8) and collect all those entries into the {v — ji) identity matrix 
lu-ix- Recall that v and /x correspond to the algebraic and geometric multiplicities of the 
zero eigenvalue of J(0), and hence there are exactly — fi) such entries to move. Such 
a movement will be accomplished by first permuting some of the first u columns in J(0) 
and then by permuting some of the first v rows of the resulting matrix. The permutations 
among the first v columns can be described by a matrix postmultiplying J(0) and we use 
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Pi to denote that matrix. On the other hand, the permutations among the first v rows 
can be described by a matrix premultiplying J(0) and we use P2 to denote that matrix. 
Thus, the matrix ^2-^^(0) Pi will be given by 

P2^(0)Pi = diag{0^,/^_^, Jn^+i(A^+i),..., JnJA^)}, 

where 0^ denotes the /x x /x zero matrix. Since Pi and P2 affect only the first v columns 

and u rows, respectively, they have the form 

Hi 1 " 

Pi^ , , (6.11) 

In—v Iji—n 

for some permutation matrices Hi and 112. 

Formally speaking, the matrix Hi describes the permutation tti given by 

TTi : (l,...,zy) {qi,...,q^), 

where 

( ni-\ hn^-i + 1, r = l,...,/x, 

Qt = < 

{ T ~ fi + a, T ^ fj. + 1, . . . ,u, 

and a e {1, ■ ■ ■ , fi} is the unique integer such that, for given r and fi, 

ni + n2-\ h Ua-i - a + j = T-ii, 

for some j G {2, . . . , 12^}- Note that, since > 1, the quantity ni + n2 + • • • + fia-i — oc 
is a nondecreasing function of a. 

Similarly, 112 is related to the permutation 7r2 given by 

7r2 : (1,...,!^) (ai,...,au), 

where 



ni H {-Ur 



a = 1, 



and p e {1, 



a — fi + p — 1, a = jj, + 1, . . . jV, 
, jj,} is the unique integer such that, for given a and fx 

ni + n2-{ h - p + s = Q!-//, 
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for some s e {2, . . . , Up}. To implement these permutatfons we let ij for j = 1, . . . ,1/ denote 
the column vectors of the standard basis in and let Hi be the xi^ permutation matrix 
whose jth column vector is Cg^ , and let 112 be the u x u permutation matrix whose kth 
row vector is ej.^ . Now observe that, if M is any v x v matrix, then the matrix 112 M Hi 
can be thought of as being obtained from M by a permutation of the columns according 
to TTi and a permutation of the rows according to tt2- 

Let us now return to the matrix F{k) defined in (6.2). By first putting it into the 
Jordan canonical form F{k) and then by applying Pi and P2 on the first v columns and 
rows of F{k)^ we form the matrix Z{k) defined as 

^A{k) B{ky 
_C{k) V{k) 

where A{k) has size /j, x n, Ti{k) has size (n — //) x (n — //), A(k) coincides with the 
submatrix of F{k) consisting of the entries in columns al and rows sn^, where a = 1, . . . , /i 
and s = 1, . . . The procedure of going from F{k) to 2{k) is similar to the procedure 
described on pp. 4638-4639 of [4], where the mappings Pi and P2 were also used. 

The small-/c limits of the block entries in the matrix Z{k) are described in the following 
theorem. 

Theorem 6.1 Assume that V in (1.1) is selfadjoint and belongs to L];(R+). Then, the 



P2F{k)Pi = P2S-^F{k)SPi, (6.12) 



asymptotics as k ^ in C+ of the matrices A{k), B{k), C{k), T>{k) appearing in (6.12) 
are given by 

A{k) ^kAi+o{k), B{k) ^ kBi+o{k), C{k) ^ kCi+o{k), V{k) ^ Vo+O{k), (6.13) 
where Ai, Bi, Ci, T>q are constant matrices, and furthermore Ai and Vq are invertible. 

PROOF: The proof for the expansions in (6.13) is similar to the proof of Proposition 4.4 
of [4] . The invertibility of T>o follows from the fact that it consists of invertible blocks and 
is given by 

Vq = diag{/,^_^, Jn^+l(An^+l), ■ • • , Jn^i^nJ}, 

where 1^-^ is the identity matrix of size {v — fx) with fj, and v denoting the respective 
geometric and algebraic multiplicities of the zero eigenvalue of J(0), and the Jn^i^a) are 
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the Jordan block matrices appearing in (6.7) corresponding to the nonzero eigenvalues for 
a = /J, + 1, . . . , K. From (6.1) and (6.10), as the (s, j)-entry of the matrix Ai we get 



where R is the matrix appearing in (6.1) and (6.2). By (5.63) and Proposition 5.7 we know 
that R acts as an invertible map from Ker [J(0)] to Ker [J(0)+]. RecaU from (6.3) and (6.5) 
that {uii,U2i, ■ ■ ■ , Ufj_j} is the basis for Ker [J(0)] and {I'lm , V2n2J • • • > '^m^mI Jordan 
basis for Ker [J(0)t]. Thus, the matrix Ai is nothing but, apart from the factor (— z), the 
matrix representation of the invertible map R with respect to the Jordan basis and the 
adjoint Jordan basis. Thus, Ai is invertible. | 

Theorem 6.2 Assume that V in (1.1) is self adjoint and belongs to L^(R+). Then, the 



asymptotics as k ^ in C+ of the inverse of the matrix Z{k) appearing in (6.12) is given 
by 



r 1 

k 



A^'[I^ + o(l)] -A^'BiV^' + 0(1) 



(6.14) 



. -V^^CiA^^ + o(l) Pq"' + 0{k) 
where Ai, Ci, and Vq are the constant matrices appearing in (6.13) and the invertibility 
of Ai and Vq is assured in Theorem 6.1. 

PROOF: The proof is exactly the same as the proof of Proposition 4.5(i) of [4]. We will 
use the decomposition formula [4,10] 



'^n—n 





'A 


B' 






C 


V 





n—jj, 



A - BV-^C 
V 



(6.15) 



Thus, as seen from (6.15) for the matrix Z{k) defined in (6.12) we have 



'A 


B' 


-1 


_C 


V 









n—jj, 



(A-BV-^C)-^ 












-BV 



-1 



'-n—ix 



or equivalently 



{A-BV-^C)-^ 



-{A-BV-^C)-^BV-^ 



-V-^C{A - BV-^C)-^ V-^C{A - BV-^C)-^BV-^ + V 
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-1 



(6.16) 



Finally, using (6.13) in (6.16) and the fact that A and T>q are invertible, we get (6.14). Let 
us note that we have written (l/k) ^^^[/^ + o(l)] in the top left block in (6.14) whereas 
that term was written as (1/A;)^^^ + o{l/k) in [4]. The two expressions are certainly 
equivalent because we can always premultiply the o(l/A;)-term by Ai^Ai- | 

We are now ready to evaluate the small-A; limit of the Jost matrix J{k) defined in 
(4.3), its inverse J{k)~^, and the scattering matrix S{k) defined in (4.6). From (6.2) and 
(6.12) we see that the Jost matrix J{k) is given by 



J{k) = /(-r , a)t [/(O, a)^]-^SP^^Z{k) P,-'S-\ 



(6.17) 



where the small- A; limit will be evaluated with the help of (3.8), (6.12), and (6.13). On the 
other hand, from (6.17) we get 



Jik)-' =SP, Z(fc)-iP25- 7(0, a)t [f{-k\ a)t]-i. 



(6.18) 



where the small-fc limit will be evaluated with the help of (3.8) and (6.14). Thus, using 
(6.17) and (6.18) in (4.6) we obtain 

S{k) = -f{k, a)t [/(O, a)t]-i5p-iZ(-fc) Z{k)-' P^S-' f{0, a)t [/(-fc, a^^K (6.19) 

In evaluating the small-A; limits of J{k), J{k)~^, and S{k), we will use a consequence of 



(3.8), namely for A; — >■ in C+ we have 

/(O, a)t [/(-r , a)t]-i = + o(l), fik, a)t[/(0, a)t]-i = /„ + o(l). 



(6.20) 



Theorem 6.3 Assume that V in (1.1) is selfadjoint and belongs to L];(R"'"). Then, as 



k ^ m C+ the Jost matrix J{k) has the behavior 

^ kAi + o{k) kBiAi + o{k) 



J{k) = sp^-' 



kCi + o{k) Vo + o(l) 



PiS-\ 



the inverse Jost matrix J{k) ^ has the behavior, as k ^ in C+, 

r 1 



J{k)-' = SPi 



k 



A^'[I^ + o(l)] -A^'BiV^' + 0(1) 
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P9S-\ 



(6.21) 



(6.22) 



and the scattering matrix S{k) defined in (4-6) is continuous at k = and we have S{k) = 
S(0) + o(l) ask^O inK with 

In 



5(0) = SP. 



-1 



'In—fi 



P2S-\ 



(6.23) 



where Ai, Bi, Ci, and Vq are the matrices appearing in (6.13), jj, is the geometric mul- 
tiplicity of the zero eigenvalue of the zero- energy Jost matrix J(0), P\ and P2 are the 
permutation operators appearing in (6.11), and S is the matrix appearing in (6.7). 

PROOF: Using (6.13) and (6.20) in (6.17) we get (6.21). Using (6.14) and (6.20) in (6.18) 
we obtain (6.22). Finally, using (6.12), (6.13), and (6.14) we obtain 

-/^ + o(l) Oik) 



Z{-k)Z{k)-^ 



-2Ci^r'+o(l) In-^. + 0{k) 



/c in R. 



(6.24) 



Then, using (6.24) and (6.20) in (6.19) we get (6.23) and S{k) = S{0) + o(l) as /c in 
R. I 



7. EXAMPLES 

In this section, we will check the validity of our formula (6.23) for some selfadjoint 
boundary conditions. 

Example 7.1 (The 5' boundary condition) As our first example, let us use the 3x3 

versions of A and B from p. S116 of [25], with A and B satisfying (1.11)-(1.13) and given 

by 





■ 1 





—a" 




-0 





-1- 


A = 


-1 


1 





, B = 








-1 




. 


-1 


. 




.0 





-1. 



where a is a real parameter and the potential V is zero. From (4.3) with /(/c, x) ~ 6**^^/3, 
we obtain J(/c), and from (4.6) we obtain the scattering matrix 5'(/c), and we have 





' —ik 


-l-\-iak- 




J{k) = 


ik 


-ik -1 


, S{k) 




. 


ik -1 . 





- i + ak 


-2i 


-2i -1 


3i + ak 


3i + ak 


3z + ak 


-2i 


i + ak 


-2i 


3i + ak 


2>i + ak 


3i + ak 


-2i 


-2i 


i-\- ak 


3i-\- ak 


3z + ak 


2ti + ak - 



(7.1) 
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From (7.1) we see that S{k) is continuous at A; = and 

12 2 



S{0) 



2 
3 

2 2 

3 ~3 

On the other hand, we see that J(0) is given by 





1 

3 



3 

2 
3 



1 
3 



(7.2) 



-1 

J(0)= 0-1 
.0 -1 

and hence its eigenvalues are Ai = 0, A2 = 0, and A3 : 



— 1, with respective eigenvectors 





-1- 




-Q- 




-1- 


till = 





, U21 = 


1 


, U31 = 


1 




.0. 




.0. 




.1. 



Thus, in the notation of Section 6 we have ni = 1, 77-2 = 1, and 77-3 = 1, with k = 3. Our 
permutation operators Pi and P2 are given by Pi = I3 and P2 = I3, and the matrix S 
appearing in (6.7) and the matrix Z[k) appearing in (6.12) are given by 



1 1- 

11, Z{k) 
.0 1. 

From (7.3), with the help of (6.13) we get 



-ik —ik {a — 2)ik 
ik —2ik —ik 
L ik -1 + ik } 



(7.3) 





—i —i 




'{a-2)i' 




Ai = 




, Bi = 


, Ci = [0 i], Vo = [-l 




i -2i 




—i 



Since fi = v = 2 and n = 3 in our example, we evaluate (6.23) and confirm that the right 
hand side in (6.23) coincides with the matrix in (7.2). 

Example 7.2 (The Kirchhoff boundary condition) The procedure in this example is 
similar to that of Example 7.1. We use the 3x3 versions of A and B from p. S117 of [25]. 
Using V = and 





-0 





1- 




--1 





0- 


A = 








1 


, B = 


1 


-1 







.0 





1. 




. 


1 


0. 
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we obtain 



J{k) = 







1 




—ik 
-1 -ik 
1 -ik 



S{k) = S{0) = 



r 1 
~3 
2 
3 
2 
L 3 



2 
3 



1 
3 
2 
3 



2 1 
3 

2 
3 
1 

3 J 



(7.4) 



On the other hand, we have J(0) given by 

--1 0" 

J(0) =1-10 

.0 10. 

The eigenvalues of J(0) are Ai = with ni — 1 and A2 = —1 with n2 = 2. We further have 



S = 



= u = 1, n = 3, Pi 

-1/V2 



-1/v^ I/V2 

.1 1/V2 



Z(fc) 



= P2 = h, 

-3ik -3ik/V2 

2V2ik -1 + 2ik 1 

. V2ik ik — 1. 



Ai^[-3i], Bi = [-3i/V2 0], Ci 







"-1 1 " 


V2i 






-1_ 



(7.5) 
(7.6) 

(7.7) 



Using the information in (7.5)-(7.7) in (6.23) we can verify that the right side in (6.23) 
coincides with S{0) given in (7.4). 

Example 7.3 (The XOR gate boundary condition) The procedure in this example 
is similar to that of Example 7.1. We use the 4x4 versions of the matrices A and B from 
[24]. Using y = and 

■i/a 
i/a i/(2o) 
i/{2a) i/{2a) 

. i/{2a) i/{2a) 
where a is a real parameter, we obtain 



B 



-0 








- 




















-1/2 


1/2 


.0 





1/2 


-1/2. 





'k/a 













-1 








0- 







k/a 





k/{2a) 







1 








J(/c) = 




{k + a)/{2a) 


, S{k)^S{0)^ 
















ik + a)/{2a) 











1 




. 





{k + a)/{2a) 


{k-a)/{2a). 




.0 





1 


0. 



(7.8) 
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On the other hand, we have J(0) given by 



-0 








- 




















-1/2 


1/2 


.0 





1/2 


-1/2. 



J(0) = 



The eigenvalues of J(0) are Ai = with m = 1, A2 = with n2 
and A4 = — 1 with 77,4 = 1. We further have 

lj, = 1/ = 3, n = A, P^ = P2 = 14^, 



1, A3 = with 77-3 = 1, 



(7.9) 



10 

10 

1-1 

1 1 J 

1/a 
1/a l/(2a) 
L 1/a 



, Z{k) 



/e/a 

k/a k/{2a) k/{2a) 

k/a 

LOO -1 J 



(7.10) 





- - 




l/(2a) 




. . 



Ai= 1/a l/(2a) , Bi = l/(2a) , Ci = [0 0], = [-!]• 

(7.11) 

Using the information in (7.9)-(7.11) in (6.23) we can verify that the right side in (6.23) 
coincides with S{0) given in (7.8). 

Example 7.4 The procedure in this example is again similar to that of Example 7.1. We 
use V = 0, and for some real parameters a and 6, we let A and B be 





-2 


1 


a" 




-0 





0- 


A = 








b 


, B = 








1 




.1 


1 


c. 




.0 





0. 



We obtain 



J{k) 



—2ik —ik —iak 
1-ibk 
—ik —ik —ick . 



, S{k) 





-i + bk ^ 
i + bk 

IJ 
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and hence 

-10 
S(0)= 0-10 
.0 IJ 

On the other hand, the matrix J(0), which given by 

-0 0- 
J(0) =001 
.0 OJ 

has eigenvalues Ai = with ni = 1 and A2 = with 77-2 = 2. We further have 

■1 0" 
1 
.0 1 0. 
' —2ik —ik —iak 
—ik —ik —ick 
. l-ihk. 

, Ci = [0 0], Vo = [l 



11 -- 


= 2, 


V 




3, 


n = 3, 


Pi- 






-1 





0- 






s 







1 





, Z{k) = 






.0 





1. 








'-2i - 


-i 






—ia 


















—i 




-i 






—ic 



(7.12) 



(7.13) 



(7.14) 



(7.15) 



Using the information in (7.13)-(7.15) in (6.23) we can verify that the right hand side in 
(6.23) coincides with S{0) given in (7.12). 
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